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Abstract 

We calculate the quadra-spectrum and quint-spectrum, corresponding to five and six 
point correlation functions of the curvature perturbation. For single field infiation with 
standard kinetic term, the quadra-spectrum and quint-spectrum are small, which are sup- 
pressed by slow roll parameters. The calculation can be generalized to multiple fields. When 
there is no entropy perturbation, the quadra-spectrum and quint-spectrum are suppressed as 
well. With the presence of entropy perturbation, the quadra-spectrum and quint-spectrum 
can get boosted. We illustrate this boost in the multi-brid infiation model. For the curva- 
ton scenario, the quadra-spectrum and quint-spectrum are also large in the small r limit. 
We also calculate representative terms of quadra-spectrum and quint-spectrum for infiation 
with generalized kinetic terms, and estimate their order of magnitude for quasi-single field 
infiation. 
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1 Introduction 



The non-Gaussian features in the cosmic microwave background radiation (CMB) and large 
scale structure (LSS) have become one of the central issues in modern cosmology. Theoreti- 
cally, a large primordial non-Gaussianity indicates new physics beyond single field slow roll 
inflation. Moreover, one can extract a number of observables from the measurement of non- 
Gaussianity, for example, the size and shape for bispectrum and trispectrum respectively. 
These observables can help us greatly in distinguishing between different inflation models. 

The three point correlation function, or the so called bispectrum has been relatively well 
investigated in the literature. Single field slow roll inflation has very small bispectrum [1]. 
Large bispectrum can be obtained from non-standard kinetic terms [2], non-local inflation 
[3], features in the Lagrangian [1], curvatons [5], old curvatons [B], multiple inflatons [7], 
quasi-single field inflation [5], multi-brid inflation [H], bifurcating trajcetores p^, thermal 
components [IT] and so on. Experimentally, the WMAP7 bound [12] for the bispectral 
estimator /tvl in local, equilateral and orthogonal shapes are given as (68% CL) 

f^f = 32 ± 21 , r^f = 26 ± 140 , /°f = -202 ± 104 . (1) 

from which one finds that a Gaussian model is still consistent with data, however the central 
values of all the three shapes are far from zero. In the near future, the Planck mission is 
able to reduce the uncertainty to A/]^^'^' ~ 5. The error bars for bispectra of other shapes 
are also expected to be greatly reduced. 

The four point correlation function, or the so called trispectrum, is also calculated recently 
in a number of models. Large trispectrum can show up in inflation with single field (for 
complete result, see [13]; for earlier work, see [lUlIS]) and multiple field [16\ inflation with 
generalized kinetic terms, curvaton scenario [T7j, multi-brid inflation [18], and so on. A (95% 
CL) bound on trispectrum from the WMAP5 data is given as [IH] 

- 3.80 X 10^ < gNL < 3.88 x 10^ , -3.2 x 10^ < tnl < 3.3 x lOl (2) 

In the Planck mission, the uncertainty can be reduced to Atnl — 560. The signal to noise 
ratio of tnl up to a given / is shown to be [20] 

(5/A^)(</)~2.2xlO-V^^/^ (3) 

Thus with a large Jnl, trispectrum is hopefully to be measured in Planck. 

For higher point correlation functions than four point, the investigations are much less in 
the literature. Some results on the consistency relations [HI [21] , and in 6N formalism [22j 
are available. However a detailed investigation is not yet present. In this paper, we shall 
calculate explicitly the five and six point correlation functions (the quadra-spectrum and 
the quint-spectrum) and extract the predictions in several models. In Section 2, we discuss 
the local shape of quadra-spectrum and quint-spectrum. We calculate the quadra-spectrum 
and quint-spectrum for single field and multiple field inflation models generically in Section 
3, for Multi-brid inflation in Section 4, and for the curvaton scenario in Section 5. We also 
calculate representative terms of equilateral quadra-spectrum and quint-spectrum in Section 
6, using single field inflation with generalized kinetic terms. We conclude and estimate the 
order of magnitude for quadra-spectrum and quint-spectrum for quasi-single field inflation 
in Section 7. 
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2 The local ansatz 



The most well studied non-Gaussian feature is the local shape non-Gaussianity. In the local 
ansatz, ({x) can be expanded locally at a spacetime point x as 

Q q 27 81 

C(x) = Ci(x) + -InlCK^) + ,^9nlCK^) + T^h^LCK^) + —^NLCK^) + 0{C!) 
5 25 125 625 

n=l 

Thus Ci's can be written in terms of the non-Gaussian estimators as 

.2 .54 ,3 . 648, .4 1944. , 

C2 — ■^/atlCi ; S3 — ; M — Y25 ' ~ T25~^^'^ ' 

The virtue of local non-Gaussianity is two folded: On the one hand, the local anzatz makes 
the analysis of non-Gaussian features simpler technically. On the other hand, in a number 
of known examples, non-Gaussianity is produced in super-Hubble scales. In this case, only 
the local shape non-Gaussianity is possible from causal considerations B 
In the Fourier space, the local ansatz corresponds to 

c(k)=ci(k)+-/^. y ^ci(ki)Ci(k-ko+-^^. y ^ci(ki)Ci(k2)Ci(k-k,-k,) 

^ 27 , f (fki f d?k2 f d^k 
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^-'I^^I^^I ^Cl(kOCl(k2)Cl(k3)Cl(k-k,-k,-k3) 

+ Oid) . (6) 
Following the standard definition, the shape functions of C(k) up to four points are defined 

as 

2 

(C(kOC(k2)) = (27r)^P(fc05'(Ekn) ' 

n=l 

3 

(C(ki)C(k2)C(k3)) = (27r)3fi(ki,k2,k3)53(5^k„) , 

n=l 

4 

(C(ki)C(k2)C(k3)C(k4)) = (27r)3T(ki, k2, k3, k,)6\J2 kn) ' (7) 



ra=l 



where -B and T denote the shape of non-Gaussianity, which are functions of the external 
momenta, as well as background quantities. Similarly, we define the shape functions for the 



•^By contract, non-Gaussianity may be produced in sub-Hubble scales, but get amplified in super-Hubble 
scales. In this case, the shape might not be local [23] . 
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five point correlation function and the six point correlation function as 

5 5 



(11 C(kn)) = (27r)3g(ki, k2, k3, k4, k,)5'iYl ' 

n=l n=l 
6 6 

(H ^C^-)) = (2^)'^(ki, k2, k3, k4, ks, k,)6%J2 ■ (9) 

n=l n=l 

In the local ansatz, the momenta dependence in the shape functions can be calculated 
explicitly, without knowing the detail mechanism how the correlation is established. Inserting 
Eq. (E]) into Eqs. ([7]), (jH]) and ©, we have 

S(ki, k2, kg) = IfNL [P{h)P{k2) + 2 perm.] . (10) 

5 

54 

r(ki, k2, kg, k4) = -gNL [P{h)P{k2)P{h) + 3 perm.] 

+TNL [P{h)P{k2)P{\k, + ksl) + 11 perm.] . (11) 



Q(ki,k2,k3,k4,k5) = ^/liVL [P(A;i)P(A;2)P(fc3)P(A;4) +4 perm.] 

+ [P{k^)P{k2)Pik,)P{\k^ + k4|) + 59 perm.] 

+ 07VL [P(A;i)P(A:2)P(|ki + k3|)P(|k2 + k4|) + 59 perm.] . (12) 

1944 

t/(ki, k2,k3, k4, ks, kg) = -Y^^NL [P{ki)P{k2)P{k^)P{k^)P{h) + 5 perm.] (13) 

+ aNL [P{k^)Pik2)Piks)P{k^)P{\ki + ksl) + 119 perm.] 

+ Pnl [P{k^)P{k2)P{k:,)P{k^)P{\k^ + k2 + ksl) + 89 perm.] 

+ INL [P{ki)P{k2)P{k^)P{\ki + k4|)P(|k2 + ksl) + 359 perm.] 

+ 5ml [P(A;i)P(fc2)P(fc3)P(|ki + k4|)P(|ki + k4 + ksl) + 359 perm.] 

+ AjvL [P(A;i)P(A:2)P(|ki + k3|)P(|k2 + k4|)P(|ki + k3 + ks|) + 359 perm.] , 



where 



324 „ , 216 
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_ 3888 _ 2916 ^ 

O^NL — -^^iNLlT-NL , PNL — "q^^NL i 

1944 ,2 , _ 1296 ,4 

INL - OnL - -^^Jnl9NL , y^NL — "^^-INL ■ U^j 

In the above calculation, we are considering only connected parts of the correlation functions. 
The disconnected parts are automatically dropped. The function Q, as well as some results 
in terms of the 6N formalism, are also obtained in [22] . 
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One should also note that Eq. (fT^ only holds when the local expansion (jlj) is valid. 
In some generalized cases (for example, multi-field inflation and multi-brid inflation, which 
we shall show below), Eq. no longer holds. Instead, the nonlinearity is expressed as a 
product of different fields at the same spatial point. In this case, the shape functions ffTOj) . 
f lTTj) . f lT2|) and f|T3l) are still the same. However, the coefficients tnl, ^nl, 4>nl, o:nl, (3nl, 
iNL, ^NL and X^L should be calculated independently. Eq. f|T^ can not be used in this 
case. 



3 Single field and multi-field inflation 

The 6N formalism ^2A\ is a powerful tool to calculate cosmic perturbations. Especially, 
it is very convenient to calculate local non-Gaussianity using the 6N formalism The 
curvature perturbation is related with the perturbation of the infiaton field as 

JV" IV'" IV"" IV""' 

C(x) = N'5<P, + —S^l + —5<Pl + -^5ct>l + -^5<Pl + 0{5<f:) , (15) 

where 50* is the Hubble crossing value of the infiaton fluctuation in the spatial flat gauge, and 
prime is derivative with respect to 0*. In the current section, we assume that nonlinearity is 
mainly developed on super-Hubble scales, which is verified in bispectrum and trispectrum. 
This assumption corresponds to assuming 

50* = (5i0* , (16) 

without non-Gaussian components. In this case, the non-Gaussian estimators can be directly 
identified with the derivatives of e-folding number as 

_ 5A^" _ 25iV"' _ l2bN"" _ l2bN""' 

^^^-^iiV^' ^^^-"^W^' '^^- 1944A^- ^^^^ 

Other estimators are given in Eq. ( |T4|) . One note from this result that as one goes to higher 
order, the results are suppressed by more powers of slow roll parameters. Here and after- 
wards, we only calculate the leading order non-Gaussianity, and neglected the contribution 
from loop diagrams, which are suppressed by at least of order 10~^, thus not likely to be 
observable. 

Note that by saying the above result is suppressed by slow roll parameters, we assume 
the time derivatives of the slow roll parameters are also small. On the other hand, if the 
time derivatives are large, for example, in the case that the infiaton field has small oscillation 
during the slow roll, the quadra-spectrum and quint-spectrum can get amplified, by powers 
of the oscillation frequency divided by the Hubble parameter. 

With the assumption that the sub-Hubble correlations can be neglected, the above calcu- 
lation is also straightforward to generalize into multi-field case, with an additional assumption 
that isocurvature perturbation is absent (we shall go beyond this assumption in the next 
section). In multiple field inflation, 5N is expanded as 
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By this expansion, we have assumed that the e-folding number difference depend only on 
the fields at horizon crossing. In other words, the uniform total energy density slice is 
the same as the uniform energy density slice for each component fields. This corresponds 
to the assumption of absence of entropy perturbation. Otherwise, one should express the 
fields on the decay hypersurface of the last field (or other hypersurface where isocurvature 
perturbation vanishes) as functions of the fields at the Hubble crossing time before doing 
the variation of 0^. The estimators f^i, Qnl, h^i and iiyL can be written as 

_ 5 E., N.jN.N, _ 25 E.,. N,,,N,N,N, 

tNL-—rT^^—^^ 9NL- ,,2n3 ' 



I'NL 



^NL = ..,,5 • (19) 



648(E.iVf)^ ' '""^ 1944 (E.AT^^' 



As Eq. f lTSj) is different from Eg. (j4]), the non-Gaussian estimators 6' atl, (pNi, C(nl, Pnl, Inl, 
6nl and Xnl are also no longer given by Eq. ( fT4l) . Instead, after taking expectation values, 
one can find 

tnl = ■ 



2^3 



e 



NL — ,^nx4 5 rNL 



(^NL = , P. 



Note that in multi-field inflation, all the relevant fields are assumed to be light. Thus the 
resulting non-Gaussianity is typically suppressed by slow roll parameters. It is interesting 
to see whether loop corrections or a hierarchy of different slow roll parameters could 
amplify the non-Gaussianities. But these aspects are beyond the scope of the current work. 
As noticed in [27], using Cauchy-Schwarz inequality 

Y,{A^A,){B,B,) > Y,iA,B,)iA,B,) , (21) 

ij ij 

where Ai and Bj are real vectors, one can derive 

TNL > ^fNL ■ (22) 

Similarly, using Cauchy-Schwarz inequality once, we have 

2Q1 fi 

Pnl > '-^9%, . (23) 



Using Cauchy-Schwarz inequality twice, we have 

^NL > ^/^. . (24) 

For other equations in Eqs. ffMj) . one can not find similar inequalities correspondingly. In- 
stead, we have checked that similar inequalities for 9^^, (pNL^ c^nl, Inl and 6nl can fail 
assuming that 6N can be an arbitrary function of 50j*. 

The quadra-spectrum and quint-spectrum should obey Maldacena's consistency relations. 
This consistency relation for general n-point correlation has already been discussed in [l4j 
and [21]. 



4 Multi-brid inflation 

As we have seen in the previous section, with the absence of entropy perturbation, non- 
Gaussianity is suppressed by slow roll parameters. However, when entropy perturbation 
is allowed, things get complicated and one might get large quadra-spectrum and quint- 
spectrum. In this section, we use multi-brid inflation to illustrate the large quadra-spectrum 
and quint-spectrum in multi-field inflation with entropy perturbation. 

We consider a simple model of multi-brid inflation: the two-brid infation with potential 

V{4>i, 02) = Vo exp (ai0i + ^202) , (25) 

Vo - li9!<Pl + gl^Dx' + 1 (x^ - x) ' • (26) 

The power spectrum and bispectrum of this model is calculated in [9], and the trispectrum 
is calculated in [18j . Here we shall generalize the calculation in [9l [18] into quadra-spectrum 
and quint-spectrum. 

In this model, inflation ends when g\(t)\ + gl'Pl, = o"^. We parameterize 0i and 02 at the 

end of inflation as ^ ^ 

0if = — COS7, 02f = — sin7. (27) 
91 92 

The e-folding number before the end of inflation is given as 



g2(^l/«l _j_ g2(^2/«2 
g20i//ai _j_ g2(/)2//a2 



+ . (28) 



where N^. is a correction of e-folding number near the end of inflation, which does not 
contribute leading order non-Gaussianity when the non-Gaussianity is large. Here 0i = (pu 
and 02 = 02* are the value at Hubble crossing. We shall omit the star subscript to avoid 
clustering of indices, as well as to keep the convention the same as that in the literature. 

To take variations with respect to 0i and 02? one need to note that 7 should also be 
written as a function of 0i and 02. The relation between 7, 0i and 02 can be established by 
the conservation of the angular field component 

01__02^01/_02/^ q'COS7 _ (TSin7 ^^g^ 
ttl 02 «! ^2 fi'l^l 92Ci2 
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Note that 0i and 02 have nearly independent and Gaussian perturbations at horizon crossing, 
which are denoted as 50i and 502 respectively. Using Eq. (^^, the e- folding number (125]) 
can be simplified into 

iV(0i,02) = (0i-0i/)/ai . (30) 
To expand Eq. (!30|) into fifth order, one have 



5N, = Ci 
6N, = C4 



gi cos 7^01 + g2 sin 7(502 
gfg2 (a250i - ai502)^ 



a^g^a^cl 
^9i92 (c+ + 5c?.) (a2(50i - ai502) " 



- ' ^ ^ 



where 



a = yal + aj , 9 = \j qI cos^ 7 + 5-! sin^ 7 , 

. fi'2a2 giai , 5'2a2 . 

c_ = sm 7 cos 7 , c+ = cos 7 H sm 7 , 

ga ga ga ga 

~ _ 92(^1 . gia2 , . 

c= sm7 COS7 . [62) 

ga ga 

One can either directly calculate the correlation functions or making use of Eq. (120 p to have 
_5c^gfgja _ 25c^c^ gfgja"^ _ ^gfgja^ 

JNL — — ^ , 9NL — 7Tr"2 A~r> ' '^NL — o fi 9 ' 

Qc+g-^a ISc^g^a^ c^g^^a^ 

I2h{c\ + hc^_)c^g\g^a^ Sc.c^gfgla^ ~^g\gW 

2164^5^3 ' - 4^V3 ' '^^^-^j^' 

_ 625c_(34 + 7c?.)c^^^^|a^ _ 3(4 + 5^)54^,6^6^4 _ gc^ 54^,6^6^4 

^Af-L — ^.o 4 fi 4 ' "^^-J- ~ 404 5 — 404 ) 



_ 3c_c35(j5(^a'^ _ 3c_c35(]'5(^a^ _ (?g\g\a^ , . 

4^;iOa4 ' 44V4 ' ^^^^"4^- (33) 

As shown in [18j, the parameter c+ can be related with the tensor to scalar ratio r and the 
spectral index ng as 

''Technically, one note that from Eq. ([29| . the zeroth order quantities satisfy 5Q!(c+ cos 7 + c_ sin 7) = 
giai, which can simplify the calculation. 
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which should satisfy c\ < 0.625 when Ug = 0.96 and r < 0.20. When r ^ 1, c+ ^ 1. 
Thus the quadra-spectrum and quint-spectrum can be large in the multi-brid model. In the 
c+ — )■ limit, the order of magnitude of the non- Gaussian estimators are of order 

iNL ~ O^NL ~ I^NL ~ iNL ~ 5nL ~ ^NL ~ C]^"^ • (35) 

One should also note that in some parameter regions, other parameters could also become 
large. 



5 The curvaton scenario 

In this section, we calculate the correlation functions up to the sixth order in the curva- 
ton scenario. The calculation is a direct generalization of [IT], in which up to four point 
correlation functions (corresponding to up to 0(f in the expansion) are calculated. In the 
calculation, we assume when the curvaton field oscillates, the potential becomes quadratic. 
We also use sudden decay approximation. One can follow [28] and [IT] to generalize the 
result to go beyond these assumptions. 

On super-Hubble scales, the curvature perturbation on uniform density slice can be writ- 
ten as 

Q{x) = 6N{x) + - / , (36) 

where subscript i can denote either infiaton (and the corresponding radiation after inflaton 
decay), curvaton or a combination of them. As there is no interaction between infiaton and 
curvaton, the infiaton perturbation (r and the curvaton perturbation are both conserved 
on super Hubble scales. 

Choosing the spatial flat slice, Eq. (!36|) for curvaton becomes 

Px = Pxe'^^ • (37) 

As usual, we assume the curvaton could have general slow rolling potential during infla- 
tion, while have quadratic potential when the curvaton oscillates. Then the initial amplitude 
X of curvaton oscillation can be written as a function of the curvaton's Hubble exit value x* 

as X = gix*)- 

At Hubble exit, the fluctuation of the curvaton field is extremely Gaussian. Thus we can 
write 

X* = X* + Six* ■ (38) 
Correspondingly, the initial amplitude of curvaton oscillation takes the form 

gix*) = 9{X* + SiX*) = 9 + y-^J , (39) 

The energy density of curvaton can be written as 

2 

(40) 
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Making use of Eq. (1371) , one can relate C,^ and 5x order by order as 



^ -'^r 6 + 12^^" a'^'"' '"'""\c^ 



>4 

9" 9" 9'" ff' " 



Cx5 = 1^24 - 60^ + 30^ + 20^ - 10^^ - 5^ + ^ J • (41) 

To relate the perturbations to we use the sudden decay approximation. Assuming the 
curvaton decays on a uniform total density hypersurface H = T, where F is the curvaton 
decay rate. Then on this hypersurface we have 

Pr + Px = P ■ (42) 

Note that on the decay hypersurface, the curvaton is oscillating, thus have equation of state 
= 0, and the radiation from inflaton has equation of state Pr = Pr/^- Applying Eq. fl5B]) . 
we have 

Pr = p.e^K-'^) , = p,e^(^^-^) . (43) 

Thus on the decay hypersurface we have 

(1 - n^)e^^'^^-'^^ + n^e3(^--^) = 1 . (44) 

We consider the standard curvaton scenario, where the primordial inflaton fluctuation can 
be ignored. In this case, ( and can be related order by order as 

^ - ^ ^2 _ 3 

Ci - <xi , ^ - ^ 

Cs _ 9, iSgY + ,,'") 9 / ^ ^ ^ , ^ ^ 3^, ^^^^ 





99"\ 


r — 


(- 


9") 


99"^ 


1 999" 


1 


9" J 


' 29" 


+ 2 



^3 ^,^2 Q^"^ ^ 

C4 27 g\?>g"' + Ag'g'" + gg"") 9 (3^?'^ + ISgg'^g" + 3^72^"2 + 4^2^,^,. 



, 9 (17^-^ + 2^^V-3^V2-4^2^V") , , / 27^7(7"\ 
+ - — + 90 + r (^-50 + ^;^J -70r - 15r , 

(46) 

C5 81 g^{lQg"g"' + bg'g"" + gg""') 135 g {6g"^g" + Qgg'^g'" + 2g^g"g"' + 9gg'g"^ + 9g^g'g"") 



Cf 16r4 ^'5 4r3 c/'S 

45 (30<7'5 + 102<7^?'3^" + 2g'g''g'" - Qg^g"g"' + 3<72^'<7"2 - 3g'g'g"") 



+ 



, 225 (3^^^ + 12^^ V + 3^ ^ 2^V(7"0 , 3 / 1000^7^7" , 45^^ (3^/.2 ^ 

+ 2r +4 ^ 1^13 + 

- 350r ( 2 + + 5r2(173 - ^^) + 630r=^ + 105r^ , (47) 

V ^ / 5 



9 



where r = 3f2^/ (4 — Using ([5]), we have 



5C2 -^Cs , _125C4 . _ 125 Cs 

6cr 54 cr 648 cr 1944 cf 



Other non-Gaussian estimators are given in Eq. (fT4l) . 

In curvaton models, large non-Gaussianity is obtained when r — )■ 0, which corresponds 
to the case that curvaton only takes up only a small fraction of total energy density when 
curvaton decays. In this case non-Gaussianity is large because curvaton has to have larger 
fluctuations in order to contribute large enough curvature perturbation in uniform total 
energy density slice. In the r — )■ limit, the non-Gaussian estimators take the form 

JNL — -L H w ' 9nl 



i-NL 



4r V g'"^ J 24r2 g'^ 

125 g'^{?,g"'^ + Ag'g'" + gg"") . 125 g^ {10g"g'" + bg'g"" + gg'' 



^NL ^ 7377-7 -r^ • (49) 



192r3 g'^ ' '''^ 384r4 g'^ 



As another special case, when the curvaton potential is always quadratic from Hubble 
exit to curvaton decay, we have g" = 0, g'" = 0, g"" = 0, g'"" = 0. In this case, the above 
result is simplified to be 



12 \r 

25/18 , 

= T7^ + 1 + 20r + 6r2 

108 \ r 

125 /54 153 

hNL = -7737777 360 + 200r + 280r2 + 60r=^ 

2592 r 

125 /675 1350 , ^ A 

iNL = + 5139 - 2800r + 3460r2 + 2520r2 + 420r^ . 50 

7776 y r J 

As is well known, in the simplest curvaton model with quadratic potential (corresponding 
to Eq. (l50ll ). in the r -C 1 limit, gNL is of order 1/r, the same scaling as /nl- Such a gNL as 
small as gNL ^ Inl is much more difficult to be measured than /nl- However, hNL scales 
as 1/r^, we have g'^1 hNL I'nl- Thus the situation for measuring hNL should be better 
than gNL, although more difficult than /nl- For similar reason, iNL is much more difficult 
to be measured than hNL- However, in the more general case , where the curvaton potential 
is not quadratic during inflation (corresponding to Eqs. (H5|l . fH6|l . (H7|) ) . there is no such 
hierarchy, and all terms are equally likely to be measured. The counter part of this effect 
for trispectrum is investigated in [25] . 



6 Non-local non-Gaussianities 

In this section, we would like to go beyond the local ansatz, which leads to the shape function 
other than those deflned in Eqs. ( fT2l) and (fT3l) . We illustrate the non-local shape using DBI 
inflaton. The complete calculation is very lengthy, with a great number of terms, which is 
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beyond the scope of the current paper. However, from the full calculation of bispectrum 
and trispectrum, we get the experience that all equilateral shape functions are similar. More 
technically speaking, the shape functions of the bispectrum generally have a large component 
of the standard equilateral shape, and have a very small mixing of orthotropic shape. Similar 
phenomenon is also observed in trispectrum. Thus we could expect that we are able to 
extract some shape information of the quadra-spectrum and the quint-spectrum, even when 
the calculation is not yet complete. 

Before a calculation, we want to estimate the order of magnitude for the result. When the 
kinetic term of infiaton is generalized, the sound speed for perturbation is generally different 
from unity. By power counting [H], one expects n-point correlation function (n > 3) scales 
as 

iC) oc P^-'/cT-' , (51) 

where is the dimensionless power spectrum, and Cg is the sound speed for perturbation. 
Thus one have 

hNL OC ~ ffjL , iNL OC ~ /^^ . (52) 

This behavior is similar with the curvaton model with non-quadratic term. 

We perform the calculation in spatial flat gauge. In the Cg <C 1 limit, the 5th and 6th 
order action of DBI inflation takes the form 

S = /'^4^[ 7a'^'/ 5aidi5(l)di5(l))6(f' ^ 3(^,(509,(50)^^0 1 
^ y \ 8c903 4cj03 8ac503 J ■ 

S = /'^4^ [ 21a''5/ 35aidMdi6(j))6(l)' ^ 15(9,^09,(50)^(5/ idMdi6(j)f \ 

In the interaction picture, the field (50 can be expanded in terms of creation and annihi- 
lation operators as 

(50^(r, k) = u{t, k)a(k) + u*{t, -k)a^(-k) , (55) 
The mode function u{t, k) satisfies the classical equation of motion 

Uk = (1 + ikcsT) e-'""''' , = ^^ctre-""'^^ . (56) 

We are fully aware that there are non-trivial transformation from Lagrangian to Hamil- 
tonian, and one also need to calculate diagrams with scalar propagation. However, to have 
some feelings of the amplitude and the shape function, we only calculate the first term each 
for (|53|) and (l5l|) . disregarding all the complexities. 

It is convent to use the commutator form of the in-in formalism 

{cm = -^ f dt'{[o{t)\Hut')]) ■ (57) 

J to 



^From the experience of bispectrum, the first term might be canceled by other terms in DBI inflation. 
However, the first term is the simplest, while hopefully carrying the general information of the shape func- 
tions. In the case of trispectrum, the corresponding term is chosen as the representive shape. Thus we still 
choose to calculate the first term. We might, with a little more effort, calculate all the terms listed in the 
action. However, an incomplete result makes no much improvement compared with a more incomplete result 
for illustration purpose. 
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For the quadra-spectrum, we have 




(C(r,ki)C(r,k2 



)C(r,k3)C(r,k4 



)C(r,k5))D 



(58) 



where = /ci + ^2 + ^3 + ^4 + ^5- For the quint-spectrum, we have 



(C(r,ki)C(r,k2 



)C(r,k3)C(r,k4 



)C(r,k5)C(r,k6))D 



(59) 



where Kq = ki + k2 + ks + k4 + k5 + kQ. Note that the quarda-spectrum and the quint-spectrum 
are proportional to and respectively, as expected. The numerical coefficients are also 
large because a large number of permutations. 

7 Conclusion and discussion 

To conclude, we discussed quadra-spectrum and quint-spectrum in inflation and curvaton 
models. In inflation models, these correlations are in general suppressed by slow roll pa- 
rameters, with the exception of oscillation, nontrivial entropy perturbation or generalized 
kinetic terms. In curvaton models, especially in curvaton models with potential other than 
quadratic terms during inflation, the quadra-spectrum and quint-spectrum can get large. 

The importance of quadra-spectrum and quint-spectrum depends on the experimental 
ability to probe them. In this paper, we simply assume the quadra-spectrum and the quint- 
spectrum are "large" when Hnl > /atl, and inl > /atl, respectively. However, this esti- 
mation need to be verifled or improved by data analysis. For CMB data, we might naively 
expect that the signal to noise ratio of quadra-spectrum and quint-spectrum are direct gen- 
eralizations of Eq. ([3]) like 



which supports our estimation when one can reach / ~thousands in Planck, with the as- 
sumption that /atl > 10. Moreover, one can hope that with the help of the future LSS 
experiments, the data points are available in three dimensions instead of two. In this case, 
the number of available data points increases, which favors higher order correlation func- 
tions. Another possibility which favors higher order correlations is that in some cases, the 
quadra-spectrum and quint-spectrum are indeed much larger than /^r^ and /^^ respectively. 
This possibility can be achieved by tuning the parameters in the models we have discussed, 
or more naturally achieved by a possibility that we shall discuss as follows. 

In quasi-single fleld inflation, which is an intermediate case between single fleld and 
multi-fleld inflation, the quadra-spectrum and quint-spectrum are large, and potentially be 
promising in probing non-Gaussianities. This is because in quasi-single field inflation, the 
non-Gaussianity in the isocurvature direction need to be projected onto inflation direction 
using transfer vertices. Using the result of |H], the non-Gaussian estimator /at^ in a turning 
trajectory model of quasi-single fleld inflation is of order 



{S/N)t 



PrfNLl' , {S/N)u ~ P^fU' , 



(60) 




(61) 
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where 9/H denotes the turning angle of inflation trajectory during one Hubble time, and 
V" denote the interaction in the isocurvature direction. To use perturbation theory in 

calculation, we need (j)/H^ ^1. As in [8j, one can estimate the order of magnitude of 
quadra-spectrum and quint-spectrum as 

hNL ~ P^'^' [e/H)' {V"'/Hf ~ {e/Hr'f%, , (62) 

iNL - P^-' {e/H)' iy"'/Hf ~ {e/HY'ft,, . (63) 

If the transfer vertex can be calculated perturbatively, then the factors {e/H)~^ and (J) /H)~^ 
are much greater than 1. Thus quadra-spectrum and quint-spectrum may play interesting 
roles in quasi-single field inflation. However, as an explicit calculation of bispectrum is 
already messy in quasi-single field inflation, one had better improve the calculation technique 
in calculating quadra-spectrum and quint-spectrum for quasi-single field inflation. 

As an extremal example, when V" ~ H ^ the isocurvature direction is completely non- 
Gaussian. In this case, if fjqi is intermediately large, say f^L = 10, we will have h^L ~ /atl, 
instead of the cubic power. 

The discussion on quasi-single field inflation can be made more general. In multi-field in- 
flation as well, non-Gaussianity might be created in isocurvature directions and be converted 
to inflation direction. Sometimes one ignores this effect not only because of simplicity, but 
also because transfer vertices are involved, which might suppress the correlation function. 
However, quadra-spectrum and quint-spectrum need less transfer vertices compared with 
/Ir^, and f^i^ respectively. Thus for quadra-spectrum and quint-spectrum for multi-field 
inflation, isocurvature perturbations also become more important. 
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